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Abstract: We study the Cramer-Rao bounds of parameter estimation and coherence 
performance for the next generation radar (NGR). In order to enhance the performance of 
NGR, the signal model of NGR with master-slave architecture based on a single pulse is 
extended to the case of pulse trains, in which multiple pulses are emitted from all sensors 
and then integrated spatially and temporally in a unique master sensor. For the MIMO 
mode of NGR where orthogonal waveforms are emitted, we derive the closed-form 
Cramer-Rao bound (CRB) for the estimates of generalized coherence parameters (GCPs), 
including the time delay differences, total phase differences and Doppler frequencies with 
respect to different sensors. For the coherent mode of NGR where the coherent waveforms 
are emitted after pre-compensation using the estimates of GCPs, we develop a performance 
bound of signal-to-noise ratio (SNR) gain for NGR based on the aforementioned CRBs, 
taking all the estimation errors into consideration. It is shown that greatly improved 
estimation accuracy and coherence performance can be obtained with pulse trains 
employed in NGR. Numerical examples demonstrate the validity of the theoretical results. 

Keywords: next generation radar (NGR); Cramer-Rao bound (CRB); Fisher information 
matrix (FIM); pulse trains; parameter estimation; coherence performance 
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1. Introduction 

Large aperture high-power phased array radar has played an important role in long-range 
surveillance, tracking and discrimination, owing to its capability of obtaining high signal-to-noise ratio 
(SNR) echoes. Typical such radars include the USA's Ground Based Radar-Prototype (GBR-P) and 
the Sea-Based X-Band (SBX) radar [1]. However, large size and heavy weight usually make them 
difficult to transport and deploy and hence, easy to be attacked in practice. In order to achieve high 
SNR gain while maintaining acceptable sensor size, a novel radar architecture has recently been 
proposed by the Lincoln Laboratory, i.e., the next generation radar (NGR) [2], where the large aperture 
phased array radar is made up of several transportable distributed sub-apertures or sub-radars. It is 
shown that NGR has improved mobility, stronger survival ability and similar processing gain 
compared with the traditional large aperture phased array radar. Moreover, an experimental NGR 
system with two radars has been constructed by the Lincoln Laboratory, which is reported in [3] to 
have obtained inspiring coherent processing gain in field tests, showing its good application prospects. 

In this paper, we consider NGR with a master-slave architecture, where all the radars transmit 
signals and only the master radar receives the echoes. It is known that the maximum echo power can 
be achieved only when we make all the transmitted signals arrive at the target at the same time and 
in-phase, namely, the coherence gain is obtained via coherent processing. However, the distributed 
architecture of NGR makes it difficult to coherently combine signals for two reasons. First, the range 
from a target to different radars may be different, leading to echoes with different propagation time 
delays and phases; second, each radar has an independent local oscillator with different transmit and 
receive (T/R) phases, which also adds phase shifts to echoes. Since both the T/R phases and the phase 
caused by propagation delay can influence the coherence gain, we add the two phases together and 
name the sum as total phase. In NGR, both the mismatches of time and phase can cause performance 
degradation. To overcome this shortage, an operation procedure of two steps has been proposed [2]. In 
the first step which is also called the MIMO mode, each radar transmits a probing signal (usually 
orthogonal waveforms) to estimate the time delay differences and total phase differences between 
sub-radars and the master radar, and they are referred to as coherence parameters (CPs) [2-5]. In the 
second step which is also called the coherent mode, all radars transmit coherent waveforms adjusted by 
the estimated CPs from MIMO mode. Clearly, the estimation accuracy of CPs greatly impacts the 
coherence gain that can be obtained by NGR, which raises two important questions: What is the best 
estimation accuracy for the CPs? How much coherence gain can we get assuming that estimation 
accuracy is achievable? 

Another problem in NGR lies in the constraint of system size, i.e., the number of radars cannot be 
arbitrarily large in practice. Thus, the maximum SNR gain that can be obtained merely through the 
spatially coherent processing of distributed radars is limited, which is unfavorable in detecting and 
tracking long-range weak targets. To settle this problem, it is natural and essential to emit pulse trains 
in NGR, which means we will accumulate the energy of echoes not only from different radars but also 
from multiple pulses. In NGR transmitting pulse trains, new questions immediately emerge: How will 
the introduction of pulse trains affect the estimation accuracy of aforementioned CPs? Are there any 
new parameters that need to be estimated? If any, what is the best estimation accuracy for those 
parameters? What is the optimal coherence performance for NGR with pulse trains? 
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A thorough review of the existing literature on NGR reveals that the present signal models in NGR 
are all based on single pulse schemes [2-5], whereas the transmission of pulse trains has not been 
considered yet. From the aspect of parameter estimation, [4] derived the CRBs of time delay 
differences and T/R phase differences for a general NGR architecture, but the CRBs of total phase 
differences are not given. Therefore, the CRBs of CPs have not been thoroughly worked out so far, 
according to the definition of CPs. In the field of performance analysis, [4] derived the performance 
bound of NGR based only on the CRBs of T/R phase differences, assuming that all time delay 
differences are ideally compensated. In [5] a formula of coherence gain taking all types of estimation 
errors into consideration was presented, but the performance bounds were not analyzed. Therefore, the 
performance bound analysis of NGR still remains an unresolved problem. 

In addition, it is worth pointing out that the parameter estimation of NGR should be distinguished 
from the parameter estimation of MIMO radar which has been studied in [6-17], despite their 
superficial similarities in emitting orthogonal waveforms. Their differences are: first, the central 
parameters of concern in NGR and MIMO radar are different. MIMO radars are focused on target 
localization accuracy, i.e., the x-y coordinate or the x-y velocity of a target, while NGR cares about 
time delay differences and total phase differences with respect to the target, i.e., the CPs. Second, the 
parameters of phases are modeled and treated differently in NGR and MIMO radar. In MIMO radar, 
phase synchronization errors are modeled as random variables which are used to evaluate the average 
performance degradation [14-16], and they need not to be estimated, thus their CRBs are of no 
interest, while in NGR the parameters of phases are modeled as deterministic unknowns that need to be 
estimated for compensation so their CRBs are of high concern. 

In this paper, we make the following contributions which also answer the questions at the end of 
paragraphs two and three. All the contributions below are useful and instructive for the system design 
and performance analysis of NGR: 

(a) The NGR signal model based on a single pulse is extended to the case of pulse trains for the first 
time, and the concept of spatial coherence is extended to joint space-time coherence for NGR. 
The extension to pulse trains benefits the detection and tracking of weak targets and helps 
control the system scale of NGR. 

(b) The original coherence parameters (CPs) of NGR are extended to the generalized coherence 
parameters (GCPs), with Doppler frequencies involved. Since target echoes coming from 
different radars usually have different Doppler frequencies, they must also be estimated and 
compensated. The extension to GCPs is essential in characterizing the multi-pulse model in (a). 

(c) The closed-form CRBs of the GCPs are derived based on the signal model in (a), and verified 
through simulations, thus providing a lower bound for the estimation accuracy of the GCPs and 
a criterion for the performance evaluation of different estimation algorithms. 

(d) The formula of coherence gain for NGR is derived and the performance bound is analyzed based 
on the CRBs in (c) with all types of estimation errors considered, thus providing an upper bound 
for the SNR gain performance of NGR. 

The paper is organized as follows: in Section 2, we present the NGR signal model with pulse 
trains and specifically define the GCPs. In Section 3, we derive the CRB for parameter estimation. In 
Section 4, we present the analytical formula of coherence performance. Simulation results and 
discussions are shown in Section 5, and Section 6 concludes the paper. 
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2. System Model and Parameter Definitions 

The system model of NGR with master-slave architecture is illustrated in Figure 1. Without loss of 
generality, we assume that there are AT radars with Radar No.l being the master radar. 

Figure 1. The master-slave architecture of NGR. 
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The pulse signal transmitted by the Ml transmitter is: 

s k (t)e j27rf ^ jd Kk = \- 



(1) 



where Sk{t) is the baseband signal of the Ml transmitter,^ is the carrier frequency, and ff k represents the 
phase of the local oscillator at the Ml transmitter. We define the effective bandwidth of Skif) as: 



ft=yl\f 2 \S t (f)\ 2 df/\\S t (f)fdf 



(2) 



where Sk(f) is the Fourier transform of Sk(i). Further, we assume that the effective bandwidths of all the 
signals are equal, i.e., fi = fik,k = \,...,K. 

As mentioned above, we consider the case of transmitting multiple pulses. Assume that TV pluses are 
transmitted consecutively and reflected by a moving point target. For convenience, assume that the 
range-rates of the target with respect to all radars remain constant and the target does not move across 
a range cell within the observation interval, which implies that no range migration occurs. The range 
from the target to the Ml transmitter when transmitting the nth pulse can be expressed as: 



R k,n ~ R k,0 ' 



■nT*R 



n = 0X"-,N-l 



(3) 



where Rk,o is the initial range from the target to the Ml transmitter, T is the pulse repetition interval 
(PRI), and R k is the range-rate of the target with respect to the Ml transmitter. 

Then, the propagation time of the ni\\ pulse in the Ml path, i.e., the path from the Ml transmitter to 
the first receiver (i.e., the master radar), can be expressed as: 



*\k,n 



Ik 
C 



0,1,--,N-1 



(4) 
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where c represents the speed of light, n^o = (^i,o+ Rk,o)/c and R lk =R { +R k represent the propagation 
time of the first pulse and the range-rate of the Ml path, respectively. 

For simplicity, we assume that the target is non-scintillating, i.e., its complex reflection coefficient 
is identical for each propagation path, and is denoted as t, . The nth pulse echo received at the master 
radar is a superposition of the returns contributed by the emitted signals from AT transmitters: 



E^cf-^y 2 *'-'"^* (5 ) 



The received signal in Equation (5) is then mixed up with the master radar's local oscillator 
e -j2Kf c t+jd x an( j (j own _ conver t e( j to base-band, where 9\ represents the phase of the first receiver. 
Substituting Equation (4) into (5), the signal after down-conversion can be expressed as: 

r un if) = t ?s k (t - V „ )e-™^A+« + Wu (, ) 



(6) 



k=\ 



where w\, n (i) is the noise of the first receiver while receiving the nth pulse. We assume that W\, n (t) is a 
temporally and spatially white zero-mean complex Gaussian random process. More specifically, 
we have: 

E{w ln (t)wlJu)} = a 2 w S(t-u)S(rn-n) (7) 

where a 2 w is a constant and d(t) is the dirac function. 

By denoting £ = g e ~ j27rf ^ +J ^ +J ^ =g r + jg. , Aff k = ff k -ff x and Az k = r^o - rn,o, Equation (6) can be 
written as: 



^h{t-T u Je J ^'" + w u {t) 



k=\ 



where Aff k represents the T/R phase difference between the Ml transmitter and the first transmitter, 
A(pl = -2jtf c ATi represents the total phase difference between the Ml path and the first path and 

(p k =27tf c TR xk lc represents the normalized Doppler frequency of the Ml path. Applying the 
previous assumption that no range migration occurs during the observation interval, we have 
Sk(t - T\k, n ) ~ Sk(t - ri£,o) approximately, and Equation (8) can be written as: 

r l Jt) = t^ k (t-T k )e JA ^ + w hn (t) 

k=l 

(9) 



= Y,^(t-T l -*T k )e jA4+M -+w hn (t) 



k=l 



where we replace n^o with ik here and hereinafter for simplicity and Ar^= ik~ t\ represents the time 
delay differences between the Ml path and the first path. 
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It is obvious that An = A# = 0 and A^f = -2nf c £±T x +A0[ = 0 , which do not need to be estimated. 
From Equation (9) we define a parameter vector consisting of deterministic unknowns as: 

v = [AT T 9 Ay T 9 y T (10) 
where Axe R (K ' l)xl , A\|/e R ( ^" 1)xl , <pe R* xl with their expressions as: 

At = [At 2 , At 3 ,-- - , At^] 7 (11) 

Ay = [A$,A#,-.,A^f (12) 

<p = [ft,p 2 ,-,?kf (13) 

Note that the time delay differences in Equation (11), total phase differences in Equation (12) and 
Doppler frequencies in Equation (13) are the GCPs as defined in Section 1. 

3. CRBs for Parameter Estimation 



In this section, we derive the CRBs of the GCPs. The Fisher information matrix (FIM) is calculated 
first. Then the CRBs are obtained by inverting FIM [18]. Finally, some remarks are given. 

3.1. FIM of Intermediate Parameters 

We find it difficult to compute the FIM of the GCPs directly, thus an intermediate parameter vector 
6 is introduced: 

6 = [^A V W,£,£] r (14) 

where t = [t 19 t 2 ,~- ,t k ] t e R Kxl . Next we compute the FIM of 6, which is denoted as 1(6), to further 

obtain the CRBs of GCPs. 

As the noise components in each pulse is independent and identically distributed (i.i.d.) Gaussian, 
the log-likelihood function of the received signal given in Equation (9) can be expressed as: 

J N-l 



\ogp(r;Q) = — L^L\r hn (t)-^s k (t-T k )e J ^-" \ 2 dt 

& w n=0 k=l 

where r = [r l 0 (t), r x x (t), • • • , r l N _ x (t)] T , and T denotes the observation time in a PRI. 
The computation of 1(6) is provided in Appendix A and the final result is: 



(15) 



1(9) = 



0 



Kx(2K+l) 



0(2^+l)x^ ^(2K+l)x(2K+l) 



where: 



A D E 
D r B F 

E r F T C 



(16) 



(17) 



(18) 
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B = ff *(*-l)(2*-l) Ig 



6 

C = OT, 



'(^-i)xi 



2 



K-l 



E = A^Tl 1 1 

^ iV |_ A (^-l)xl A (^-l)xl J 



1 ~ L A ^xl A ^xlJ 



5. 



(19) 
(20) 
(21) 

(22) 
(23) 
(24) 



where 1^ represents the K identity matrix, 0 p x q and l p x q represent ap * q zero matrix and ap x q 
matrix with all elements being unit, respectively. We denote l p x q as J p when p = q. Note that T, A, B 
and C contain the second-order derivatives with respect to t, A\|/, (p and [£ r , &] T , respectively, D 
contains the second-order derivatives with respect to both A\|/ and (p, E contains the second-order 

derivatives with respect to both A\|/ and [£ r , &] T , and F contains the second-order derivatives with 
respect to both (pand [£ r , £J r . For convenience, we define the CRB matrix (CRBM) as the inverse 

matrix of FIM and denote CRBM by CRB e , i.e., CRB e = I _1 (8). 

3.2. CRB of Time Delay Differences 

The CRBM of t can be obtained from Equations (16) and (17): 

CRB T = T 1 

I, 



8n 2 j3 2 N\{\ 
1 



2 K 



(25) 



$n 2 /3 2 NSNR in K 

where SNR in =\£\ 2 / a 2 w is defined as the input SNR at the receiver. 

Next we apply the chain rule of CRB to compute the CRBM of Ax by using the relationship 
between Ax and x: 



which gives a Jacobian matrix as: 



A ^=^-^i ? k = 2, — ,K 



3|W r i j 1 

3 T L (K-l)xl *K-\_\ 



(26) 



(27) 



where g(x) = At. Then the CRBM of At is: 
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CRB. =Ml) CRB , ^ g(T) ^ 



1 



"■(K-l)xl 



"-K-1 



8n 2 j3 2 NSNR ir 
Sx 2 j3 2 NSNR in 

Extracting the diagonal elements of CRBa t , the CRBs of time delay differences are expressed as: 

1 1 1 



(28) 



2 n2 ' 



N SNR in 4x 2 p 



k = 2, — ,K 



(29) 



which coincides with the result in [4] for the case of K transmitters and one receiver, 
except that a coefficient of l/N is multiplied, indicating the improvement on estimation accuracy by 
pulse integration. 

3.3. CRB of Total Phase Differences and Doppler Frequencies 
For clarity, we rewrite matrix G as: 



G _ 2 



s u 
V c 



where: 



A D 
D r B 



U = 

v = u r 



(30) 

(31) 

(32) 
(33) 



Using the matrix inversion lemma [19], the inverse matrix of G, which corresponds to the CRBM 
of A\|/, (p and [£ r , £] r , can be expressed as: 



G 



2 



S U 

V c 

(E-UC'V)" 1 -IT'l^C-VST'lTf 1 

-crVcs-uc-V) -1 (c-ve-'u)- 1 



(34) 



where the submatrix (E - UC V) 1 corresponds to the CRBM of A\|/ and (p which are of interest and 
(C - VL -1 U) _1 corresponds to the CRBM of &f which are nuisance. 
From Equations (21) and (31-33), we have: 
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E-UC V = 



A D 
D r B 



A' IT" 
V D' 



KN 



[E T F r ] 



(35) 



where: 



A' = A — EE r =|^| 2 



^K-l ^ K-l 



1 T IH2 N(N-l) 

U' = D EF = \£\ — i '- 

KN 11 2 



-ll I -lj 



1 T . M2 N(N-l)(2N-i) 
D = B FF = \a — ^ ^ 



KN 



* 2(2N-\)K K 



V = V' T 



(36) 
(37) 

(38) 
(39) 



Then, using the matrix inversion lemma once again, we have: 



(E-UCT'V) 



A' U' 
V D' 



(A' - U'D'" 1 V')" 1 - A'-'U'CD' - V'A'-'U')" 1 
-D'-'V'tA'-U'D'-'vy 1 (D'-VA'-'U')" 1 

where (A - U'D'^V')" 1 and (D' - V'A'^U') -1 correspond to the CRBM of Ay and (p, respectively. 



(40) 



Using the following two equations: 



(I r +cJ r ) '=I r — J r 

V K K) K l + cK K 



\ ( 



„ ^ K-l J K-l ^-K-l „ J K-l ^-K-l T r ^ K-l ^ K -I 



K 



K 



J V 



K 



(41) 
(42) 



1 



where c is a constant and c * , we obtain: 

K 



1 (2N-1), 

CRB. = -^ Ul K x +3 Kl ) 

Av SNR in N(N+1) k ~ x K ~ lJ 



CRB„ = 



1 



9 SNR in N(N 2 -l) K 



(43) 
(44) 



where the definition of SNR in is the same as in Equation (25). 

Finally, the CRBs of total phase differences and Doppler frequencies are expressed as: 

1 2(2N-l) 

CRB ^= -, r. k = 2,-,K 

SNR m N(N+l) 



(45) 



and: 
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1 6 



SNR in N(N 2 -l) 



(46) 



respectively. 

3.4. Remarks on CRB Results 

Based on Equations (29), (45) and (46), we have the following remarks on the CRB results: 

(1) All the CRBs are irrelevant to the number of radars K, which is the characteristic of the master- 
slave architecture. In the estimation procedure of MIMO mode, K orthogonal signals are 
extracted from the mixture echo of the master radar to estimate K—\ time delay differences, K-\ 
phase differences and K Doppler frequencies. This implies that more radars bring more 
parameters to be estimated, and the estimation accuracy does not increase as K increases. 

(2) The CRB of the total phase differences in Equation (45) only depends on input SNR and the 
pulse number N. Note that the CRB of the T/R phases in [4] is proportional to the squared ratio 
of the carrier frequency to the effective bandwidth ifjp) . Obviously, the latter is much higher 
than the former under the assumption of narrowband signals, and performance analysis based 
merely on the mismatch of T/R phases would be inappropriate and more or less discouraging, as 
we will see in Section 5.3. 

(3) When the pulse number TV is relatively large, the CRB of Ar^, A(p{ and ^is proportional to 1/N, 
4/N and 6/N , respectively. Intuitively, for the estimation of Aifc the coherent integration of TV 
pulses is equivalent to increasing the input SNR TV times. In contrast, for the estimation of A(p{, 
due to its coupling with the estimation of cpk, the coherent integration of N pulses results in a 
reduced equivalent SNR gain of N/4. The CRB of cpk descends the fastest among all the GCPs as 
TV increases. 

4. Coherence Performance Analysis 

In this section, we analyze the coherence performance of NGR with pulse trains. The signal model 
in coherent mode is presented first. Then the formula of coherence gain is provided. Finally, some 
remarks are given. 

4.1. Signal Model in Coherent Mode 

We assume that the GCPs are stable when NGR switches from MIMO mode to coherent mode, so 
that the estimates for GCPs in MIMO mode can be applied to adjust the phases and time delays on 
transmit. The estimates for GCPs are defined as: 



AT = [A£ 2 ,A£ 3 ,---,A£^f 



(47) 



Av = [A&A&-.,A&f 



(48) 



(49) 
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We use 5t k - Ar h - Af k , 8(j)[ - A0* k - Afy and 8(p k =q> k -q> k to represent the estimation errors. In 
coherent mode, the nth pulse transmitted by the kth transmitter can be expressed as: 

^(^ + AfJ^' 2 ^ +7 ' A ^" 7 ' A ^"^' w (50) 

where s(t) is the baseband coherent waveform, s(t + Ai k ) represents the time delay adjustment, e~ JA(pk 
compensates the total phase difference and e~ J ^ k ' n compensates the phase caused by Doppler. 

From Equations (8), (9) and (50), the nth pulse echo received by the master radar after propagation 
and down-conversion is: 

£=1 

= 2 J ^(r-r 1 -A^+ArJe l ; +^(0 (51) 



Since all the GCPs have been compensated on transmit, the N pulses received by the master radar 
can be directly accumulated to achieve coherent integration, and the final integrated signal is: 

N-l 

K0 = 2>i,«W 



n=0 

N-l K N-l 



*s(t-T l -ST k )e j ^"'+Y J w hn (t) (52) 

n=0 k=l n=0 

= tf i s(t-T l - Sr k )e m A{8cp k ) + £ w Kn (t) 



k=l n=0 

N-l 



where A{dq> k ) = ^ e jS(Pk ' n represents the impact of Doppler estimation errors. 

n=0 

For simplicity, dik, Scpl and Scpk are assumed to be independent with Gaussian distribution, i.e., 
dz k ~ N(0, c 2 ), k = 2,..., K, dfy - N(0 9 a 2 0 ), k = 2,-, K, dcp k ~ N(0, a 2 ,), k = 1,.., K. Note that 5x x = dfr = 0 
and the mean square errors (MSE), i.e., a 2 , q\ and a 2 are lower bounded by the CRBs derived in 
Section 3. 

4.2. Coherence Performance Analysis 

It is obvious from Equation (52) that all the three estimation errors, i.e., dtk, Sfy and Scpk, will 
degrade the coherence gain. For convenience, we assume that linear frequency modulation (LFM) 
signal with a large time-bandwidth product is adopted and the range envelope after pulse compression 
is modeled as a sine function: 

p(t) = smc(7TBt), \t\<T p (53) 

where B is the bandwidth of LFM with /? 2 = 5 2 /12 and T p represents the pulse width. 

First we calculate the averaged power of the output signal. The signal in Equation (52) is sampled at 
t = x\ where a peak after integration is expected, and the noise-free sampled signal can be expressed as: 
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r(T l ) = &p(ST k )e J ^A(S<p k ) 



k=l 



(54) 



where the sine function p(i) is applied. 

The computation of the averaged power in Equation (54) is provided in Appendix B and the final 
result is: 



so 



= ^[|r(r 1 )| 2 



Hfl 2 



[i+(*-i)r 2 ]. 



N-\ 



N + 2^(N-n)*e 



+ 



2(K-\)T x .e 2 +(K-l)(K-2)T l 2 »e~' 7 * 



n=0 



-|2 



(55) 



where T^o 2 ) = £[p(<Jr 2 )] and r 2 (o*) = £[[p(<Sr 2 )| 2 ]. Note that T^o 2 ) and r 2 (o 2 ) are both monotone 
decreasing functions of a? where ri(0) = ^(0) = 1 and ri(oc) = 72(00) = 0. When a 2 is relatively small, 
e.g., 0 < g 2 t < 0.03, T\ and T 2 can be approximated by Taylor expansion of p(t): 

pit) =sin(xBt)/(xBt) 



3! 



5! 



l{jcBf) 



(56) 



6 120 



However, when a 2 is large, the lower-order Taylor expansion is no longer suitable and polynomial 
fitting is used to calculate T\ and r 2 . In detail, the curves of T\(a 2 ) and r 2 (o\) are obtained by Monte- 
Carlo simulations first and then fitted into two groups of polynomial coefficients. In summary, the 
analytic formulas of Ti(o 2 T ) and T 2 (o 2 ) are expressed as: 

' — J 0<cr r 2 <0.03 

(57) 



polyval[P x , B 1 o) ) , 0.03< a) < 40 



and: 



l + 2c,.<7 2 , 



0<ct 2 <0.03 



polyval[P 2 , B 1 g) ) , 0.03< a] < 40 



(58) 



where c\ = -{kB) 16, c 2 = (ttB) /40 and polyval(P, x) represents the value of a polynomial evaluated atx 
with vector P containing the polynomial coefficients. In the simulations, we use the following two 
groups of 20-order polynomial coefficients: 

5.8018e-24 -2.4016e-21 4.6149e-19 -5.4646e-17 4.4633e-15 -2.6672e-13~ 

1.2070e-ll -4.2245e-10 1.1577e-08 -2.4994e-07 4.2574e-06 -5.7041e-05 

5.9680e-04 -4.8194e-03 2.9538e-02 -1.3430e-01 4.3949e-01 -9.9565e-01 

1.4904e+00 -1.4167e+00 9.8609e-01 



(59) 
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P 2 = 



1.4338e-23 -5.8985e-21 1.1259e-18 -1.3233e-16 1.0719e-14 -6.3475e-13 

2.8435e-ll -9.8387e-10 2.6616e-08 -5.6629e-07 9.4847e-06 -1.2461e-04 

1.2739e-03 -1.0004e-02 5.9224e-02 -2.5756e-01 7.9401e-01 -1.6530e+00 

2.1792e+00 -1.7018e+00 9.2949e-01 



(60) 



where the polynomial coefficients are stacked in row-order and descending powers. 

From Equation (52), the noise power after integration is Pmo = Na 2 w , then the formula of coherence 
gain for NGR with pulse trains can be expressed as: 



'NGR MSlave 



_SNR out _P so /P t 
SNR in P SI /P f 



NO 



]_ 

~N 



N-l 



N + 2^(N-nye 



2 9 



n=\ 



n=0 



(61) 



[i+(*-i)r 2 ] 

+ 2(K-l)T^e 2 +(K-\)(K-2)T?.e~ cj; 

where P SI /P M = \^ j C 2 W represents the input SNR. 
4.3. Remarks on Coherence Performance 



Based on Equation (61), we have the following remarks on the coherence performance: 

(1) If a? = q\ = q\ =0, i.e., the estimators are ideally accurate, the maximum SNR gain of ^ 2 7Vcan 
be obtained. Note that K is the ideal SNR gain of a master-slave single-pulse NGR consisting 
of K radars, and TV is the ideal SNR gain of integrating TV pulses coherently for a master-slave 
NGR with a single radar. This concept of joint space-time coherence indicates that more pulses 
can be employed to exchange for fewer radars to obtain a desired coherence gain, which makes 
the NGR system more flexible. 

(2) If a 2 = q\ = q\ = oc, i.e., the estimation accuracy is extremely low, the minimum SNR gain of 1 
or 0 dB can be obtained. The explanations are as follows: a 2 = oc means that the echoes emitted 
from other radars can be hardly aligned with the echo emitted from the master radar, and thus no 
spatial coherence gain can be obtained. Meanwhile, a 2 9 = oc means that a random Doppler 
compensation phase is multiplied to each transmitted pulse of the master radar, so the coherency 
of the N pulses is completely corrupted and the signals are integrated incoherently instead of 
coherently, which amplifies the power of signal by a factor of N. Note that the noise power is 
also amplified by TV times, we have the minimum SNR gain of 1. 

(3) If we replace a 2 , a 2 * and a 2 9 with the CRBs given in Equations (29), (45) and (46), respectively, 
then (61) gives an upper bound for the coherence gain of NGR versus input SNR. 



5. Numerical Results 



In this section, numerical results are presented and discussions are conducted to verify the CRBs of 
the GCPs and evaluate the coherence gain formula. 
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5.1. Simulations on CRB 



Since the CRBs are not impacted by the number of radars K, for simplicity we consider a NGR 
system with master-slave architecture consisting of two radars. The LFM signal with a time-bandwidth 
product of 100 is applied in simulation and all the GCPs are estimated by the maximum likelihood 
estimator (MLE), which is asymptomatically unbiased and efficient, using Monte Carlo simulations 
with 500 iterations per SNR value. 

The MSE and CRB of the time delay difference versus SNR are shown in Figure 2(a), where N = 2 
and the bandwidth B increases from 10 MHz to 100 MHz. It is seen that larger B gives better 
estimation accuracy. As SNR increases, the MSE descends quickly and approaches the CRB in the 
high-SNR region. Since the time delay measurement in radar system is usually scaled by the temporal 
resolution, i.e., 1/5, we measure the CRB of the normalized time delay differences which is, from 
Equation (29), given by: 

3 1 __L_ v i v 



CRB n 



Note that the CRB B A Tk in Equation (62) is irrelevant to /? and B. Figure 2(b) depicts the MSE and 
CRB of the normalized time delay difference for B = 100 MHz with different number of pulses. We 
see that the MSE under each TV value asymptomatically approaches the corresponding CRB curve, and 
the CRB is decreased by 3 dB when TV is doubled. 

Figure 2. The logarithmic MSE and CRB of (a) the time delay difference and (b) the 
normalized time delay difference vs. input SNR. 
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The MSE and CRB of the total phase difference and the first Doppler frequency (we have two 
Doppler frequencies to estimate since K = 2) with different number of pulses are shown in Figure 3(a) 
and (b), respectively. As expected, the MSEs approach the CRBs for high SNR in all cases. The CRB 
of Doppler frequency declines faster than that of total phase difference when TV doubles, as we have 
analyzed in Section 3.4. 
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Figure 3. The logarithmic MSE and CRB of (a) the total phase difference and (b) the 
Doppler frequency vs. input SNR. 
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5.2. Simulations on Coherence Performance 

The upper bound of SNR gain in Equation (61) vs. input SNR is plotted in Figure 4 with different 
number of pulses when K=2. For a fixed N, the SNR gain ascends from 0 dB to a maximum limit value 
as the input SNR increases. When TV doubles, the maximum gain increases by 3 dB, i.e., more pulses 
means better performance. Moreover, the input SNR required to achieve a desired SNR gain can be 
reduced by applying more pulses if we have a fixed number of radars, as can be seen from Figure 4. 
Finally, it should be pointed out that the SNR gain curves plotted here only provide a performance bound 
for NGR, since the practical estimation algorithms cannot reach the CRB in the low- SNR region. 

Figure 4. SNR gain vs. input SNR with different number of pulses. 




0 5 10 
SNR jp (dB) 



5.3. Comparative Simulations 

As mentioned in the second remark of Section 3.4, the CRB of the total phase differences in 
Equation (45) is irrelevant to /? and f C9 while the CRB of the T/R phase differences derived in [4] is 
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proportional to (fjfi) . The CRBs of these types of phase errors and the SNR gain based on each of 
them are shown in Figure 5(a) and (b), respectively, where N = l,f c = 1 GHz and the bandwidth B 
increases from 10 MHz to 500 MHz. It is seen that when B is relatively small, e.g., 
5=10 MHz, the CRB of the T/R phase differences is nearly 50 dB higher than that of total phase 
differences and an input SNR of about 60 dB is required to reach the optimal SNR gain of 6 dB, which 
makes NGR seem impractical with narrowband signal. As mentioned in Section 3.4, performance 
bound based merely on the T/R phase errors, as done in [4] is inappropriate, for the misalignment of 
time delay also plays an important role in the performance analysis of NGR. Total phase differences 
which consist of both the T/R phases and the phase caused by time delay, characterize the signal model 
more accurately. 



Figure 5. Comparisons between total phase differences and T/R phase differences: (a) the 
logarithmic CRB (b) SNR gain vs. input SNR. 
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6. Conclusions 



We have extended the NGR model based on a single pulse to the case of pulse trains, so that the 
coherence gain can be obtained from both spatial and temporal integration. Accordingly, the coherence 
parameters (CPs) in [2-5] are extended to the generalized coherence parameters (GCPs), with Doppler 
frequencies involved and the total phase differences introduced. Based on the signal model of NGR 
with master-slave architecture in MIMO mode, the closed-form CRBs of the GCPs are derived and 
verified using simulations. Moreover, in order to investigate the impact of estimation errors of GCPs 
on coherence performance, we developed the signal model of NGR in the coherent mode and derived 
an analytical bound of SNR gain, which relies on the CRBs of all the GCPs. Our simulations show that 
the introduction of multiple pulses in NGR not only improves the estimation accuracy and the 
maximum coherence performance, but also reduces the input SNR required to achieve a desired SNR 
gain, compared with the case of a single pulse. 
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Appendix A 

In this Appendix, we develop the FIM for the intermediate parameter vector 6 defined in 
Equation (14), based on the log-likelihood function given in Equation (15). It is easy to get the first 
derivatives of the GCPs as (The primitive function log p(r; 6) and the integration time T are omitted 
for simplicity): 



d 2 p 

= — -Re 



d 2 p 

3M < 



Re 



d 2 p 
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n=0 V k=l 
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(63) 

(64) 

(65) 
(66) 

(67) 



Applying the properties of orthogonal waveforms, we have the second derivatives of the GCPs as: 
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(68) 
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From Equations (68) to (77), the FIM defined in Equation (16) is obtained. 
Appendix B 

In this Appendix, we develop the averaged power of the output signal. From Equation (54), we have: 
= E[\r(T l )\ 2 ] 



1 SO 



= E[r(?y(jS\ 

[\_k=l J L/=l 

= Iff E\j^p{dT k )p{8T^- m A{d^M\^Pi) 



(78) 



k=\ 1=1 



The i.i.d. property of the estimation errors and the following equations will be applied hereinafter: 



St, = Stf = 0, E[A(Sq> k )] = E[A\Sq> k )], E[e ±jS * 
the sum in Equation (78) in the following three cases: 



= e" 2<T ' and E[e ±m ~\ = e~ 2<7 " . Then we calculate 
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This part of the sum can be expressed as: 
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This part of the sum can be expressed as: 
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In addition, we have: 
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From Equations (79) to (83), Equation (78) can be simplified as: 
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